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The precision of double-beta (/3/3) decay experimental half-lives and their uncertainties is reeval¬ 
uated. A complementary analysis of the decay uncertainties indicates deficiencies due to small size 
of statistical samples, and incomplete collection of experimental information. Further experimental 
and theoretical efforts would lead toward more precise values of /3/3-decay half-lives and nuclear 
matrix elements. 


PACS numbers: 23.40.-s,21.10.-k,28.20.Ka,29.87.+g 

Double-beta decay was proposed by M. Goeppert- 
Mayer p] as a nuclear disintegration with simultaneous 
emission of two electrons and two neutrinos 

( Z , A ) —y (Z + 2, A ) + 2e + 2z/ e . (1) 

There are several double-beta decay processes: 2/3“, 
2/3 + , e /3 + , 2e and possible decay modes: two-neutrino 
(2i/), neutrinoless (Ozz) and Majoron emission (y°) 

(Z, A) (Z ± 2, A) + (2e ± ) + (2P e , 2i/ e or X °). (2) 

The /3/3-process has been extensively investigated in 
the last 80 years m is]. These efforts have led to obser¬ 
vations of the two-neutrino decay mode and deduction 
of decay half-lives. It is the rarest presently-observed 
nuclear decay. In a recent analysis of /3/3-decay data, 
A.S. Bar abash [4 claimed that we could deduce precise 
values of experimental half-lives and extract the corre¬ 
sponding nuclear matrix elements (NME). Unfortunately, 
these claims are premature and not very beneficial for 
the field. The main goal of this work is to investigate the 
present status of /3/3-decay research, reanalyze the data, 
and produce realistic assessments and recommendations 
that could accelerate overall progress. 

The /3/3-decay T^ 2 are available from multiple sources 
0|5]. A brief summary of the recent National Nu¬ 
clear Data Center (NNDC) evaluated or adopted values 
is shown in Table [I] and plotted in the lower part of Fig. 
[]] This plot contains 12 adopted half-lives for nuclei of 
practical interest. The NNDC evaluation is completely 
based on standard U.S. Nuclear Data Program proce¬ 
dures, and its validity has been extensively scrutinized 
using theoretical predictions and nuclear systematics ar¬ 
guments [5|. 

To gain a complementary insight on the statistical 
sample size, as shown in the Table [l| I will resort to 
a non-traditional approach and consider Benford’s Law 
[6]. Figure [I] illustrates that relatively large statistical 
samples for thermal neutron capture cross sections and 
B(E2)t adopted values of 421 and 427 nuclei [7, 8], re¬ 
spectively, are in agreement with the law. The Pear¬ 
son’s cumulative test statistic values for nuclear reaction 
and structure data sets are 8.584 and 15.0, respectively. 


TABLE I: Adopted /3/3-decay T 2 / 2 for 0 + -4 0 + transitions. 
Data are taken from the Ref. [5]- 


Parent nuclide 

2?/a(y) 

A v /U +X \y) 


(4.39±0.58)xl0 19 


76 Ge 

(1.43±0.53)xl0 21 


82 Se 

(9.19±0.76)xl0 19 


96 Zr 

(2.16±0.26)xl0 19 


100 Mo 

(6.98±0.44)xl0 18 


116 Cd 

(2.89±0.25)xl0 19 


128 Te 


(3.49±1.99)xl0 24 

130 Te 

(7.14±1.04)xl0 20 


136 Xe 

(2.34±0.13)xl0 21 


130 Ba 


(1.40±0.80)xl0 21 

150 Nd 

(8.37±0.45)xl0 18 


238 u 


(2.00±0.60)xl0 21 


These cumulative values could be compared with upper- 
tail critical values of y 2 distribution [9] with 8 degrees 
of freedom. A similar analysis for the adopted values of 
/3/3-decay T fj 2 0 is not possible because one is supposed 
to have at least 5 counts in each of the nine bins flO] . 
or the number of /3/3-decay observations has to quadru¬ 
ple. Consequently, there is no sense in applying it to the 
presently-available /3/3-decay half-lives due to the limited 
number of observed transitions and large experimental 
uncertainties. 

From the dawn of the /3/3-decay era researchers knew 
the importance of complete experiments when both re¬ 
leased energy and angular distribution of decay products 
have been recorded. The first direct observation of two- 
neutrino decay mode in 82 Se m and subsequent obser¬ 
vations in 150 Nd and 48 Ca [T2l[T3] were made using such 
techniques. These discoveries employed the time pro¬ 
jection chambers (TPC) that contained small amounts 
of target material and, consequently, generated limited 
statistics. The experiments were very complex; however, 
observation of two-electron events has provided clear ev¬ 
idence of the double-beta decay process. 

Further developments have progressed using advanced 
commercially-available detectors [HH3]. and large quan¬ 
tities of enriched isotopes. Unfortunately, the usage of 
commercial detectors often leads to incomplete exper- 
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FIG. 1: Benford’s distributions for adopted values of thermal 
neutron capture cross sections, B(E2)t and /3/3-decay T^y 2 
data sets [3 El E] 


iments when only energy release information has been 
collected. To illustrate this point, we should consider /3/3- 
decay search in 76 Ge. The lower part of Fig. [2] depicts 
a chronological record of the measurements [14]. The 
different half-life values with large uncertainties for two- 
neutrino mode of decay are often explained by the rela¬ 
tively high background in the earlier experiments. These 
experiments were based on an essentially-single source of 
enriched 76 Ge isotope and slightly different detector fab¬ 
rication and shielding technologies. All of these measure¬ 
ments relied heavily on an excellent energy resolution of 
Ge detectors and suffered from the lack of electron track¬ 
ing information. Consequently, a final /3/3-decay spec¬ 
trum could contain the contribution of single-electron 
events. High energy resolution is absolutely essential for 
observations of neutrinoless of /3/3-decay, when decay will 
produce a sharp peak corresponding to a Q-value between 
parent and daughter nuclei. Unfortunately, the quantum 
world is very diverse, and background processes may af¬ 
fect the results. It has been demonstrated recently that 
7 -ray transitions in Pb and 76 Ge could produce a 2039 
keV signal and obscure the decay signature nsuis]. 

Electron tracking information is crucial in /3/3-decay re¬ 
search. An observation of two-electron events in addition 
to energy release information would help to suppress the 
radioactive background and decisively prove /3/3-decay 
process in a particular nucleus. Perhaps, in addition to 
HPGe technologies, it could be interesting to remeasure 
76 Ge T ^ 2 in a NEMO experiment that is based on a 
more extensive set of observables [ 20 ], or explore germa¬ 
nium tetrafluoride (GeF^. Similar scenarios have been 
developing with /3/3-decay search in 136 Xe [2TH23] . 

In light of this disclosure, it becomes clear that excel- 
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FIG. 2: Chronological record of three different groups of mea¬ 
surements: Hubble constant US], T 1/2 of 6 He(/3) D3 and 
76 Ge 2 "03/3) Hi]- These graphs are courtesy of the cited pub¬ 
lications. 


lent energy resolution of HPGe detectors alone may not 
be sufficient for an observation of neutrinoless mode in 
76 Ge. Therefore, the preset day state-of-the-art GERDA 
and Major ana experiments mug may not be able to 
provide an ultimate answer due to the incomplete nature 
of the data recording. 
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Analysis of the lower part of Fig. [ 2 ] shows that the 
“precise” T ^ 2 value of Barabash [4] strongly contradicts 
the latest result of the GERDA experiment [14]. This 
example shows that it is too early, at this point, to state 
a high precision of adopted half-lives and NME, and large 
error bars of the previous NNDC evaluation m are more 
appropriate. All NNDC values were produced from the 
experimental half-lives when new measurements would 
trigger the data reevaluation every 5-7 years. 

A complementary analysis of the upper and lower parts 
of Fig. [2] shows that the presently-discussed situation 
with /3/3-decay measurements is not unique in physics. It 
is rather a common occurrence when initial, pioneering 
measurements are not very accurate and often discrepant. 
It definitely has happened to Hubble constant and 6 He 
/3-decay experiments mm, as shown in the upper and 
middle part of Fig. [2] This figure content is based on 
original graphs borrowed from the Refs. 

in m hi]. 

Other recent cases of inconsistent data include the dis¬ 
crepant decay scheme and half-life of 139 Ba [25] . and cross 
section values. In the recent review of neutron cross sec¬ 
tion deficiencies, M.B. Chadwick m compiles an im¬ 
pressive list that includes 235 U and 197 Au fast neutron 
capture cross sections. The last cross section value was 
used in calibration of the stellar nucleosynthesis KaDO- 
NiS database [27] and has a broad impact across neutron 
physics. We are surrounded by a large number of im¬ 
perfect nuclear data sets and constantly work on their 
improvement. 

The experimental /3/3(2z/)-decay half-lives include con¬ 
tributions from the nuclear structure effects and decay 
kinematics. T 2 J 2 values are often described as follows 

T^O+L 0 + ) = G^.Z)|Mr T -|M|"P, (3) 

where the function G 2u (E,Z) results from lepton phase 
space integration and contains all the relevant constants 
[1|. Equation ^ highlights a direct dependency of ex¬ 
perimental NME on the calculated values of phase space 
factors. 

Both Barabash’s and NNDC’s evaluated NME were 
based on the best available phase space factor calcula¬ 
tions in 2010 and 2013, respectively [2, 29. [30]. Table 
|TT| shows the evolving values of 76 Ge phase space fac¬ 
tors, and the table data raise a question about calcula¬ 
tion limitations and possible model dependency of the 
PSFs. The observed discrepancies between recent cal¬ 
culations of Kotila & Iachello [30] and Stoica & Mirea 
m create reason for concern. The exact values of phase 
space factors are needed in order to deduce the precise 
values of experimental NME for comparison with recent 
theoretical calculations of Senkov and Horoi [32] . It will 
be highly beneficial for the field if a third group of theo¬ 
rists will clarify the situation. The formula <© gets even 
more complex for the neutrinoless mode where half-life 


TABLE II: 76 Ge /3/3-decay phase space factor (PSF) values 
in 10 -21 yr _1 for 0 + —> 0 + transition. 


Authors 

Year 

PSF 

Boem V Vogel [2] 

1987 

130.54 

Doi et al. [28] 

1993 

53.8 

Suhonen V Civitarese [29] 

1998 

52.6 

Kotila V Iachello |30| 

2012 

48.17 

Stoica V Mirea [31] 

2013 

43.9 


also depends on a neutrino mass, and neither neutrino 
mass or NME can be easily disentangled. 

Finally, the double-beta decay experimental half-lives 
have been reanalyzed using standard nuclear data tech¬ 
niques. The analysis of the two-neutrino mode of /3/3- 
decay data sets indicates their uncertainties are high due 
to relatively small experimental statistics and incomplete 
collection of decay information. We live in an era of /3/3- 
decay observations with large uncertainties. It may take 
years before future complete or “almost” complete ex¬ 
periments like low-pressure TPC mm or simultaneous 
measurement of positrons and 7 -rays [33] will improve 
the data. The precision of experimental NME strongly 
depends on the quality of phase space factors, and addi¬ 
tional theoretical calculations are necessary to clarify the 
discrepant results. 

Double-beta has been an exciting nuclear physics phe¬ 
nomenon that has slowly revealed its experimental signa¬ 
tures in the last 30 years. It may take many painstaking 
experimental and theoretical efforts before the process 
will be well understood and measured. It is not an un¬ 
usual situation in a history of science; in fact, it is a 
very common case. Therefore, one can assume that it 
may take another 30 years before all decay modes, ex¬ 
perimental half-lives and NME values will be finalized. 
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